Proceedings of the VIII Serbian-Bulgarian Astronomical Conference (VIII SBAC) 
Leskovac, Serbia, May 8-12, 2012, 

Editors: M. K. Tsvetkov, M. S. Dimitrijevic, K. Tsvetkova, O. Kounchev, Z. Mijajlovic 

Publ. Astron. Soc. "Rudjer Boskovic", 2012, Belgrade 



MODELING OF MOLECULAR CLOUDS 
WITH FORMATION OF PRESTELLAR CORES 

SAVA DONKOV 1 , ORLIN STANCHEV 2 , TODOR VELTCHEV 2 , 

o i 

On \ 1 Department of Applied Physics, Technical University 

£^ ■ E-mail: savadd@tu-sofia.bg 

^ | 2 Department of Astronomy, Faculty of Physics, University of Sofia 

E-mail: o_stanchev@phys.uni-sofia.bg, eirene@phys.uni-sofia.bg 

. 

Abstract. We develop a statistical approach for description of dense structures (cores) in 
molecular clouds that might be progenitors of stars. Our basic assumptions arc a core mass- 
' density relationship and a power-law density distribution of these objects as testified by nu- 

merical simulations and observations. The core mass function (CMF) was derived and its slope 
in the high-mass regime was obtained analytically. Comparisons with observational CMFs in 
several Galactic clouds are briefly presented. 
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1. INTRODUCTION 



Molecular clouds (MCs) are the typical regions of star formation in galaxies. Re- 
cent high-resolution observational studies in the Milky Way reveal that MCs exhibit 
an extremely complex, clumpy and often filamentary structure (e.g. Andre et al. 
2010, Mensh'chikov et al. 2010), with column and spatial densities varying by many 
orders of magnitude. The detected large non-thermal linewidths which scale with the 
size of the cloud or of its larger substructures (e.g. Larson 1981, Solomon et al. 1987, 
Bolatto et al. 2008) have been interpreted as indicators of the presence of supersonic 
turbulence. Numerous works in the last two decades have demonstrated that this 
supersonic turbulence is among the primary physical agents regulating the birth of 
stars. It creates a complex network of interacting shocks, where dense cores form at 
the stagnation points of convergent flows. Thus, although at large scales turbulence 
can support MCs against contraction, at small scales it can provoke local collapse 
of the emerging prestellar cores. Hence, the timescale and efficiency of a protostar 
'^_J , formation depend strongly on the wavelength and strength of turbulent driving source 

. C?.' (Klcsscn, Heitsch & Mac Low 2000, Krumholz & McKee 2005). 

An important structural parameter in analytical and semi-analytical models of 
star formation in MCs (e.g. Padoan & Nordlund 2002, Hcnnebelle & Chabricr 2009, 
Vcltchcv, Klcssen & Clark 2011) is the probability density function (p-PDF), which 
gives the probability to measure a given density p in a cloud volume dV. As demon- 
strated from many numerical simulations, its shape is approximately lognormal in 
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isothermal, turbulent media that are not significantly affected by the self-gravity (e.g. 
Vazqucz-Scmadcni 1 994, Padoan, Nordlund & Jones 1997, Ostrikcr, Gammic & Stone 
1999, Federrath, Klessen & Schmidt 2008). Its lognormality should correspond to the 
same feature of the observed probability distributions of the column density N (N- 
PDFs) in MCs, due to the correlation between the local values of p along a single line 
of sight (Vazquez-Semadeni & Garcia, 2001). 

On the other hand, it has been argued that the PDF displays scale-dependent 
features and/or its shape evolves significantly in time (Federrath, Klessen & Schmidt 
2008, Pineda et al. 2010). The lognormality is typical in the low-density, predomi- 
nantly turbulent regime, whereas at higher column densities a power-law tail is emerg- 
ing. Such high-density power-law (PL) tail is a characteristic feature of iV-PDFs in 
evolved MCs where star formation processes already occur or just start (Kainulainen 
et al. 2009, Froebrich & Rowlcs 2010). That is confirmed as well from analysis 
of numerical simulations of clouds dominated by gravity (Ballesteros-Paredes et al. 
2011). A consistent theory of cloud structure must take into account the character- 
istics sketched above. Probing the N- and p-PDFs can be used to set up constraints 
to analytical star formation theories. 

In this work a statistical approach is suggested to derive the mass function of high- 
density (prestellar) cores that are possible progenitors of stars. Our starting point is 
a description of the PL tail of the p-PDF that is representative for the MC regions 
populated by these objects. 

2. STATISTICAL DESCRIPTION OF DENSE CORES 

2. 1. DESCRIPTION OF THE HIGH-DENSITY POWER-LAW TAIL OF THE PDF 

Kritsuk, Norman & Wagner (2011) showed that the p-PDF evolves from purely log- 
normal shape to a combination of a lognormal 'hat' and a PL tail. Schematic rep- 
resentation of this distribution is given in Fig. 1. We use a standard designation 
of the logarithmic normalized density: s = \g(p/po), where pq is a normalization 
unit. Parameters of the PL tail are the lower (si ow = lg(pi/po)) and the upper limit 
(s up = lg(p2/po)) of the high-density range and its slope q. Thus, the PL tail is 
described by: 



where the coefficient A s — (q In 10)/(10 9Sup — 10 9Slow ) is obtained from the requirement 




(1) 



f£ dP > = 1 - 
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2. 2. CORE MASS-DENSITY RELATIONSHIP 

In observational studies, the internal structure of MCs is usually described by a sam- 
ple of discrete condensations, delineated through different procedures and labelled 
'clumps' or 'cores'. In our approach, the set of individual clumps in a considered vol- 
ume is represented by an ensemble of statistical, not individual objects as the physical 
characteristics of both groups obey the same statistical relations (e.g. size-mass rela- 
tion, density distribution). Then the p-PDF corresponds to the density distribution 
of statistical clumps. Those with densities in the range of the PL tail are dense cores 
and hereafter we label them just "cores" for simplicity. 

Our basic physical assumption is the existence of a core mass-density relationship: 



lg(^)=,lgP) (2) 

where mo is an arbitrary unit of normalization and the power-law index x is negative 
and assumed to be fixed within the whole PL tail (Lombardi, Alves & Lada, 201). 
Further we adopt the natural presupposition about a statistical relation between core 
masses to, densities p and sizes I, 



m po V l 



and obtain by use of Eq. [2] core size-density and size-mass (with size normalization 
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\TOn/ 1 — X \/n/ 



unit Zo) as well: 

Is/ 

,p J 1-x 

(-)^lg(f 

\mo/ 1 — x \to< 

Taking into account the one-to-one correspondence between the core quantities p, 
m and I, one derives from Eqs. [TJ [H and 0] the probability distributions of core masses 
and sizes: 

dP(m) = A m (— ) qX d\g(— ) , (6) 

dP(0 = 4(rV**te(^ . (7) 



where A m and A; are calculated from the normalization of probability measures. Note 
as well the obvious relation: dP s = dP(m) = dP(l). It defines the distributions of 
core masses and sizes in the PL tail which we label hereafter m-PDF and /-PDF, 
respectively. 

It seems natural to choose the mean density of the cloud (p) and a fixed fraction 
K of its size L as normalization units (Veltchev, Klessen & Clark, 2011): 

Pa = (p) , Io = kL (8) 



2. 3. AVERAGED CORE QUANTITIES 

The average core density could be calculated in two alternative ways: 
Arithmetic average: 

F) -^/(-)-(-) 9 ^(-) = (^r)(-) 

Pi 

Logarithmic average: it is a generalization of the geometric average in case of con- 
tinuous density distribution and is defined as: 



l>2 




Analogically, one can define and calculate arithmetic and logarithmic average of core 
mass ans size by use of Eqs. [6] and [7] It is important to note that the approach of 
logarithmic averaging leads to a physically natural relationship: 



(P2/ P l) q+1 - 1 

{P2/P1Y -i 



(9) 



(HL) =(£.) .(L) 

Vm / lg V / lg ^*o' 
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The transformation from dimensionless to physical quantities is independent on the 
averaging approach: p = po (p/po), Tn = mo (m/mo), I — lo Q/lo)- Recalling Eq. [3j 
one obtains: 

mo m 



m 
~pl? 



po io p r 

We point out that the second sign of equality in this equation holds only in case of 
logarithmic averaging (cf. Eq. [TTj) . 



(12) 



2. 4. RELATION BETWEEN THE NORMALIZATION UNITS AND THE FORMULA FOR 
TOTAL NUMBER OF CORES 

The core mass is a fundamental non-observable quantity and it is crucial to avoid ar- 
bitrariness of choice of its normalization unit. Therefore we derive a relation between 
po, lo and mo using the requirements for volume and mass conservation considering 
the whole PL tail. Its total volume Vtot is calculated as follows: 



-PDF 



i 3 m = i 3 n N tot 



l-PDF 



O-'Uot 



I \ 3 

T ) dP(i) , 



where TV; is the number of cores with size I and N tQt is a measure of the total number 
of cores. The limits l\ and li correspond to the density limits p\ and p2 according to 
Eq. 2] and one obtains after integration: 



Yu 



qx 
1-x 



+ 1 



W/Oo) 



q+- 



(Pi/Po) 



q+- 



{P2/P0Y - {Pl/P0) q 



= r Q N tot -Q{q,x) . (13) 



Analogically, the total mass of the cores is calculated through: 



M tot = ^2 mN m = ... = m N tl 



qx 



i-PDF 



qx + 1 



(P2/P0) q+ * ~ (Pl/P0) q " 
(P2/P0) q - (Pl/P0) q 



(14) 



In view of M to t = P • Kot = Po ■ {p/ Po) • Vtot, w e get the relation between the 
normalization units: 



Po 



E° (£. 

Pl \P0 



qx + 1 



lk + 1 



(pi/pi) 



9+ J 



1 



(P2fpi) q+ i 



1 



(15) 



Referring to Eq. ITBIand the relation: rpoV = rM = Mtot = Po • {pl Po) ■ Vtot , where 
M and V are the mass and volume of the whole cloud, we obtain the formula for total 
core number as follows: 



JVt„t = -v 



^PoJ 



(16) 



where the filling factor r accounts for the mass fraction of cores within the whole 
cloud and the density average could be arithmetic as well logarithmic. 
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3. DERIVATION OF THE CORE MASS FUNCTION 

3. 1. TOTAL NUMBER OF CORES AND SCALES DEFINED BY CHOSEN DENSITY 
THRESHOLD 

We consider an MC as a hierarchical set of spatial scales. They are defined as effective 
sizes of subregions (or a set of subregions) in the cloud delineated through chosen 
density thresholds. Thus, a fixed density p' corresponds to a scale which contains 
the mass of the whole PL-tail substructure with densities in the range p' < p < pi. 
Since the mass-density power-law index x is negative, the large the density of a core 
the lower its mass. To derive the cumulative CMF correctly, one must take into 
account the contribution of each spatial scale in the hierarchical structure defined 
by threshold p' . The latter is a product of the total number of scales N sca i es (p') 
contained in the volume over the threshold p' and the total number of cores N to t{p') 
at each scale. Repeating the procedure in the previous sections for a subsection of the 
PL tail [p',p2], one obtains N tot (p') according to Eg. [TBI while r, (p/po) and Q(q,x) 
are functions of the given p' (instead of p\): 



Ntot(p') 



r(p') 



— j(p') ■ Q(q,x,p') 

Pq> 



(17) 



Letting p 1 <C p 2l it follows in both approaches of averaging that (p/po)(p') oc (p'/po) 
(Eqs.[51and [TUl) and Q(q, x, p') oc (p' / po) < ' 1 ~ x ^ x (Eq.[J31 in this context it is important 
to mention that simulations predict q < —1.5). In this case N to t(p') oc (Kp'V^ 3 ) ' 
(p'/po)~ 1/x . 

The total number of scales N scll \ cs (p') is easily calculated from the filling factor 
r and the requirement of mass conservation at each density threshold p' . The total 
mass of cores over the given threshold p' in MC with mass M and volume V is 
po ■ (p/po){p') ■ Vtot(p') = M tot (p') = r{p')M = r{p') Po V whence 



Vtot (p') = r(p') 



V 



Obviously rp V = M tot = N sca \ cs (p')-M t ot(p') = N sca i es (p')-p -(p/po)(p / )-V to t{p / ) 
and after simple transformations one gets finally: 



jy; 



scales 



(p') = 



rV 



i<i> 



(pO • Vtot(p') 



r(p') 



(18) 
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3. 2. THE CUMULATIVE CMF 



The cumulative CMF M(m') is derived by counting all cores with masses in the range 
m' > in > m2, corresponding to densities p' < p < p2- In our case, one has to 
multiply the total numbers of cores and scales over threshold p' = po(m' /mo) x (Eqs. 
QUandCEBl): 



Af(m') = 



£(") 

K 6 \p J 



i(-) (19) 



(P' < P2) 



Note that the slope Y of a power-law differential CMF is exactly the same of its 
corresponding cumulative function. In other words, we derived a differential CMF 
with slope —1, typical for fractal structures (Elmegreen 1997). 

In view of the high densities in the PL tail, we may assume that most (if not 
all) cores are gravitationally unstable and contract in timescales given by the free-fall 
time Tff oc p~~ x / 2 . Then a time- weighted CMF would be more representative for results 
one can expect from observations. Such time weighting can be done introducing a 

weighting factor of each scale proportional to Tg 1 oc (p'/po) oc (m'/mo) 1 / 2 . Then 
the time- weighted cumulative CMF will have a slope modified by addend x/2: 



= £(-) 

k \mo/ 



1 \ -l+l 



(20) 



(p' < Pa) 



4. DISCUSSION 

Most recent studies dedicated to the CMF show that its high-mass slope is close 
or identical to that of the initial stellar mass function Y = —1.35 (Salpeter 1955). 
In Table [1] we illustrate the predictive capability of our model in comparison with 
CMFs, derived in some observational works and from simulations. The derived slopes 
and their variations could be explained by variations of the core mass-density power- 
law index in the range > x > — 0.7 (with a single exception with an extreme value) 
which is consistent with the typical values of this quantity for the inner parts of several 
MCs as derived by Donkov, Vcltchev & Klessen (2011) and testified by comparison 
with the observational study of Lombardi, Alves & Lada (2010). 

These preliminary results are stimulating to develop further our model, including 
results from recent numerical simulations and/or analytical estimates of the time 
evolution of the p-PDF (Girichidis et al. 2012). 

Acknowledgement: T.V. acknowledges partial support by the Deutsche Forschungs- 
gemeinschaft (DFG) under grant KL 1358/15-1. 
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Tabic 1: Slopes of CMFs derived from various authors in comparison with our model 
predictions. 



Galactic MC 



Rcf. 



Note 



Slope of the CMF 



Pipe 1 

Orion 2 

Orion A 3 

Perseus 4 

Ophiuchus 5 

Perseus, Serpens, 6 
Ophiuchus 

Simulations 7 
(PP, PPV) 



(a) lognormal 

(b) time-weighted 

time- weighted 



-1.35 

-1.35 

-1.3 ±0.3 

-1. ±0.1 
-2.15 ±0.08 

-1.35 

-1.3 ±0.4 



-0.7 

-0.7 

-0.6 ±-0.6 

~ 
- -2.3 

-0.7 

-0.6 ±0.8 



-1.15 < T < -1.35 -0.3 > x > -0.7 



[1] Alves, Lombardi & Lada 2007, [2] Nutter & Ward- Thompson 2007, [3] Ikeda & Kitamura 
2009, [4] Curtis & Richer 2010, [5] Andre et al. 2007, [6] Enoch et al. 2008, [7] Smith, Clark 



& Bonnell 2008 
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